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ABSTRACT 

Sequences of infinitesimally thin, uniformly rotating, self-gravitating relativistic discs 
with internal two-dimensional pressure have been constructed. It is shown that in 
weaker relativistic configurations the sequences undergo a continuous bifurcation from 
a disc to a ring structure, while in stronger relativistic cases the sequences terminate 
at the mass-shed limit where gravitational forces are exactly balanced by centrifugal 
forces. 
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1 INTRODUCTION 

Relativistic, self-gravitating discs may have played an im¬ 
portant role during the process of galaxy formation in the 
early universe. For example, they can be considered as pos¬ 
sible predecessors of galaxies today that contain a massive 
black hole in their centre. Also, during the final evolution of 
a close binary system, consisting of two neutron stars, the or¬ 
bit shrinks rapidly due to the emission of gravitational waves 
and the two stars merge finally into one object. Shortly after 
this coalescence, the system can be approximately described 
as a highly flattened disc-like structure which may display 
relativistic effects. These systems may then be modeled, in a 
simplified approach, by infinitesimally thin discs where the 
vertical variation has been averaged onto the midplane. 

Uniformly rotating relativistic discs with zero internal 
pressure have been studied in detail by Bardeen & Wagoner 
(1971). Here an expansion method was used in which the 
problem was reduced to a set of ordinary differential equa¬ 
tions, which were then solved numerically for the higher or¬ 
der corrections. Nevertheless, using this method, strongly 
relativistic disc structures could also be calculated accu¬ 
rately. Recently, the analytical solution of this problem in 
terms of integral equations was obtained (Neugebauer & 
Meinel 1993), which could subsequently be solved utilizing 
hyper-elliptic functions (Neugebauer & Meinel 1995). This 
solution represents the first analytic solution of a continu¬ 
ously extended rotating object in general relativity. In these 
approaches a given equation of state (p = 0) and rotation 
law (0 =const.) are assumed, and the Einstein-equations 
are solved directly. 

The method of mirror images is another frequently used 
option to construct models of relativistic discs. Here one be¬ 
gins with a known axisymmetric metric (e.g. Kerr), makes 
a cut in the upper half of the p — z plane where p and 2 are 
cylindrical coordinates, and identifies the metric functions 


there with those of a symmetric cut in the lower half which 
joins the first at the 2 = 0 plane (see Pichon & Lynden- 
Bell 1996, and references therein). This method yields disc 
solutions where the corresponding matter properties, i.e. 
the energy-momentum surface distribution, must be inferred 
from the pasted solution after the two cuts have been iden¬ 
tified. The method was first used for relativistic discs by 
Turakulov (1990). 

It is known that, in the Newtonian limit, the uniformly 
rotating pressure-less dust disc is unstable to all perturbing 
modes (Binney & Tremaine 1987), a result which most likely 
carries over into the relativistic regime (Bardeen & Wagoner 
1971). The addition of pressure tends to stabilize the disc. 
For Newtonian discs it was shown recently (Kley 1996) that 
uniformly rotating, pressurized discs undergo a continuous 
bifurcation process into a ring-like structure via a sequence 
of dumb-bell shaped equilibria. This process has its three- 
dimensional analogue in the bifurcation of rotating stars into 
rings (Eriguchi & Sugimoto 1981). 

In this paper we generalize the relativistic zero-pressure 
disc solutions by imposing a given (polytropic) equation of 
state relating the surface pressure with the surface density. 
This pressure may be thought of as being generated by the 
random internal peculiar velocities in a galaxy, and it has 
a stabilizing effect. This work on the other hand also gen¬ 
eralizes the Newtonian results on thin Maclaurin discs and 
their bifurcations to the relativistic case. 

In Section 2, the relevant equations are given and the 
numerical method is outlined in Section 3. In Section 4 the 
results of the numerical computations are presented and in 
Section 5 we conclude. 


2 BASIC EQUATIONS 
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2.1 Einstein equations 


We shall treat the problem of an infinitesimally thin gaseous 
disc with internal pressure acting only in the plane of the 
disc. The problem is treated in full general relativity. The 
disc is assumed to rotate rigidly with the angular velocity 
17. We use a cylindrical coordinate system 

X\ = P, X 2 = Z, X3 = X4 = t ( 1 ) 

where the disc is located in the equatorial plane (z = 0). 
The line element is written in the form used by Bardeen & 
Wagoner (1971) 

ds 2 = e 2/1 ( dp ' 2 + dz 2 ) + p 2 B 2 e~ 2u (dip — ujdt) 2 — e 2v dt 2 . (2) 

This metric describes general axisymmetric and stationary 
configurations in general relativity. Four metric functions 
(p, iz, B, and w), which depend only on the coordinates p 
and z, have to be determined. 

For the energy-momentum tensor we use the expression 
for an ideal gas. We assume relativistic units and set G = 
c = 1 in all formulae. 

Tab - (e + p)u a ub + p {gab - . ( 3 ) 

with 

c)z 

Na = ~dx°- = ( °’ lj °’ 0) ' (4) 

In equation (|ii]), e is the total mass-energy density and p 
is the gas pressure. These are three-dimensional quantities. 
The last term, containing a projection along the z-axis, has 
been introduced to ensure that the pressure acts only in the 
equatorial plane, i.e. the forces perpendicular to the disc (in 
the z-equation) vanish. 

The four velocity in case of pure rotation (no radial or 
vertical motions) is given by 

u = m 4 (0, 0, n, 1) 

where 17 is the angular velocity as seen by an observer at 
infinity. Throughout this work we assume rigid rotation, and 
thus 17 is a constant. From the normalization of the four 
velocity (u a u a = — 1) one obtains with the given metric 


4 e v 

(i-^) i/2 ’ 

(5) 

where 


v v = pPe _2l/ ( 17 — u>) 

(6) 


is the rotational velocity in the locally non-rotating frame 
of reference. 

The Einstein equations G ab = 47r T a b yield four second 
order partial differential equations for the four metric po¬ 
tentials (see also Bardeen & Wagoner 1971) 


-2fi 


-V(BViz) - —e~ 4v p 2 B 2 X7u> ■ Viz 
B 2 


47T 


(e+p) 


1 + Vcp 

1 — v 2 


+ P 


(7) 


e ~2(a+/j) 

2 P B 2 


[V(p 2 P 3 Vw) - 4p 2 B 3 Vis ■ Vw] = 
-8n(e + p)v lp (l - u 2 ) \ 


( 8 ) 


—?rV(pVB) = 8 np, (9) 

pB 

and 

- [-4pB(p,, pp +p,zz) +4Ve‘ / • V (Bpe~ v )+ 

e~ 4v p 3 B 3 Vu> ■ Vw] = 87re. (10) 

Here the V-operator is used in cylindrical coordinates (p, z). 

Since the disc is located in the z = 0 plane, it can 
be treated mathematically as a boundary condition to the 
otherwise empty space. The above equations reduce in the 
case of a vacuum to 


V(Vz/) 

= ip 2 P 2 e“ 4l 'Vw ■ Vw - IVB • X7v 

2 B 

(11) 

V(p 2 Vw) 

2 

=-Va • Vw, 

a 

(12) 

V(pVP) 

= 0 

(13) 

A c p 

= ip 2 P 2 e _4 "Vw ■ Vw + 4vP ■ Viz 

4 B 



1 dv 



- Viz-V v+~ 
p dp 

(14) 


where A c denotes the Cartesian Laplace operator (for the 
coordinates p, z), and 


B 3 



In the disc we define the two-dimensional total surface 
energy density a, and the surface pressure P through 

e = a(p)S(z)e- 2 ", p = P{p)8{z)e~ 2 ^ (15) 

where <S(z) denotes the usual Dirac delta function. 

The proper density and pressure in the disc are then 
given by 

<T P = cre -,i , P p = Pe -M . (16) 

From momentum and energy conservation T ab ; b = 0 we 
obtain the radial hydrostatic equation in the plane of the 
disc 


P p ,p + (u p + Pp)U' p — 0, (17) 

where U' is defined through 
e 2U ' = ( u 4 ) 2 » (l - viy 1 e~ 2v . 

In deriving ( 0 ) we assumed a constant rotation rate 17. 


J = 2ty 


The total angular momentum of the disc is given by 

V: 


(v + p) 


ip — 2u t ~)2 2 i 

-—e B p dp 


1 — V 2 


and the total gravitational mass by 
M = 2n (a + 2P)Bpdp T 217J 

J Disc 


(18) 


(19) 


2.2 Equation of state 

To close the set of equations, thermodynamic relations have 
to be added. We present calculations for two different equa¬ 
tions of state. In the first, more explorative, test case we 
assume a simple relation between the total proper surface 
energy density er p and the proper (surface) pressure P p 
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(20) p, 2 = —2na. 


(29) 


where 7 denotes a polytropic exponent which is set to 3 for 
all models. This first equation of state has been used only 
for the non-rotating models. 

For the non-rotating and all rotating models we as¬ 
sumed an equation of state that obeys, in analogy to three- 
dimensional configurations, an isentropic relation for the 
proper two-dimensional surface quantities 

dcr p - ap + 0 Pp dap = 0, (21) 

° P 

where cr® denotes the proper surface rest mass density. The 
the total surface density cr p is given by 
0 . i 

CTp (Jp + &pi 

where <jp denotes the proper (surface) internal energy den¬ 
sity. In both cases equation ( |17| ) can be integrated and yields 
the condition of hydrostatic equilibrium. For the isentropic 
case one obtains 

e" -\/l — v% ap \ = const. (22) 

Up 

For the isentropic equation of state we assume further a 
polytropic relation 

Pp = K (up) 7 = (7 — l)up. ( 23 ) 

In the non-relativistic limit the metric potential v turns 
into the Newtonian potential and the hydrostatic equation 
reduces to the well known relation 

1 7 PI 2 0 2 , 

v H- r 12 = const. 

7-1 a 2 

For 7 = 3 there exists an analytic solution in the Newtonian 
limit, the Maclaurin disc solution (cf. Binney & Tremaine 
1987). I11 this case the rotation rate f2 and the polytropic 
constant K are related through 



where A' max is the maximum value of K which is taken for 
zero rotation f2 = 0, and f2 c is the maximum rotation rate 
for the pressure-less dust disc which has K = 0. The surface 
density is given by 

<r(p) = &a\Jl- (p/pd) 2 , (25) 

where op denotes the central surface density of the disc and 
Pd is the disc radius. 


2.3 Boundary conditions 


In the equatorial plane (2 = 0) the potentials are continuous. 
The presence of the disc however is given by jump conditions 
for the first derivatives of the metric potentials. By vertically 
integrating equations (|^to^) the following jump conditions 
are obtained: 


V,z 


U,z 

Ell 

B 


2-7T (<7 4- P) 

—8iv(o + P) 


1 + 1$ 
1 — i$ 

fl — U! 

1 — 1$ 


+ P 


4nP 


(26) 

(27) 

(28) 


These equations form the boundary conditions within the 
disc region p < pd- Outside, for p > pd all vertical derivatives 
with respect to 2 vanish. On the polar axis (p = 0, 2 > 0) 
the radial derivatives of all potentials are zero. Because of 
local flatness the additional relation 

B = e" +p (30) 

holds on the axis, and will be used in obtaining the numerical 
solution. At infinity the potentials v, ui, p vanish, and 5 = 1. 


3 METHOD OF SOLUTION 

The vacuum Einstein equations ( fi~i] - [i~i| ). together with the 
boundary conditions ( p6| - |i^) in the disc, form a very com¬ 
plicated set of non-linear partial differential equations and 
we resort to a numerical solution. As it is numerically incon¬ 
venient to work with boundary conditions at infinity, we first 
perform a coordinate transformation that takes infinity to 
a finite distance, say 1. Specifically, we transform the equa¬ 
tions into the new coordinates u and v which are defined 
by: 

u = —-—, and v= —-—. (31) 

P + Pd z + Pd 

The disc has the radius Ud = 0.5 in these coordinates, which 
is held fixed in all models. 

To obtain the actual numerical solution, the partial dif¬ 
ferential equations (llr|l4j) are discretized on a suitable grid 
of NU x NV grid-cells. In the models presented here, a res¬ 
olution of 128 x 128 grid-cells was used; only in the stronger 
relativistic cases was a higher resolution of 256 x 256 used. 
The radial grid-points are usually spaced equidistantly in 
physical space within the disc region p < pd (u < 0.5) 
and equidistantly in u outside this region. The vertical grid- 
points are spaced equidistantly in v. Specifically, the first 
radial grid-point u 1 lies at u = 0 (p = 0) and the last one 
unu+i at u = 1 (p = 00), and similarly for the u-coordinate. 
The discretization leads to a set of four coupled systems of 
linear equations which have a banded structure. The non¬ 
linear right hand sides of eqs. ( fll[|F^ ) are not linearised, but 
rather updated at each cycle during the process of iteratively 
solving the matrix-equations. This iterative solution of the 
matrix equations is obtained by the successive over relax¬ 
ation (SOR) method, which is easy to program and allows 
a simple treatment of the boundary conditions. 

The boundary conditions at infinity and on the axis 
are implemented straightforwardly. At infinity the values of 
v, ui, and p are set to zero and B is set to 1. On the axis and 
outside the disc the derivatives vanish which is implemented, 
for example, to second order in the grid spacing as 

vij = (3 v 3 j - P2j)/4, 

where the index j denotes the vertical grid-points. Only the 
actual jump conditions within the disc (v = 0, u < ltd) have 
to be treated with more care. Starting from the minimum 
radius of the disc (ring) and then moving subsequently to 
larger values, the updated surface density is calculated from 
the hydrostatic equation (^) using the old values of the 
metric potentials. The jump conditions j26| - |29| ) are then 
used to obtain the new values of the potentials in the disc. 
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Figure 1. The angular velocity ft and the total gravitational 
mass M as function of central redshift z for the pressure-less dust 
disc. The symbols refer to the numerical solution and the lines to 
the analytical solution. 


Using this procedure, the boundary values are updated at 
each iteration cycle of the solver of the matrix equation. 

The procedure of the solution depends on the type of 
solution sought. For the cases of disc and ring, and also 
for the intermediate ‘dumb-bell’ case which connects these 
two, slightly different methods have to be applied (see Kley 
1996). In all cases the outer radius pd of the configuration is 
held fixed at the value Ma = 0.5 in dimensionless units. Ad¬ 
ditionally, the central potential v c = v{ 0,0) is fixed during 
the iterations which is the equivalent of fixing the central 
redshift z, which is given by 


\/—gtt(0, 0) 


(32) 


where 


gtt = (pBe v uj) 2 - e v . 


(33) 


The problem is then to determine the rotational angular ve¬ 
locity ft and the pressure constant K simultaneously during 
the iterations. In the disc regime, either ft or K can be fixed, 
and the other parameter can then be calculated from the hy¬ 
drostatic equation (p2|). In the case of fixed O, the constant 
in (^) is fixed at the disc radius pd and K is then obtained 
from evaluating ( p^ ) at the centre of the disc. In regions 
where K(fi) may be double valued, which occurs in fact in 
the vicinity of the bifurcation to a dumb-bell sequence, it 
may be advantageous to switch the procedure and fix K. In 
this case fl is calculated from (p2|) at a point within the disc, 
as at the origin the rotational velocity is zero. In the disc re¬ 
gion both methods lead to the same results, which provides 
a good test for the internal consistency of the method. The 
latter method of holding K fixed during the iterations can 
in fact used up to that point at which the central density of 
the disc has reached zero. This dumb-bell region can also be 
calculated by fixing the central density <r(0) and obtaining 



p/p« 

Figure 2. Comparison of numerical and analytical results for 
three metric functions (see text) for the dust disc having a central 
redshift of z = 2.97. Lines refer to the analytical solution and 
symbols to the numerical solution. The radial coordinate p is given 
in units of the disc radius p<j. 


fl and K simultaneously from the hydrostatic equation Jiil ) 
at two different locations. Again, this is consistent with the 
previous method. The dumb-bell part ends with a model 
which has zero central density at p = 0, and from here on 
the ring solutions continue. In this sequence, in addition to 
the outer radius pd, the inner radius p\ of the ring is also 
held fixed. As before, K and fl are obtained from the hy¬ 
drostatics. The global accuracy of the method can be tested 
by comparing the asymptotic behaviour of v and c o with 
the total mass and angular momentum of the disc. As the 
spherical radial coordinate r = (p 2 + z 2 ) 1 ' 2 —> oo 


M 



(34) 

(35) 


Numerically, the evaluation of M and J at large distances 
agrees with the values obtained by integration over the disc 
© and to within one percent. Only for stronger rel¬ 
ativistic discs the agreement becomes weaker, due to the 
strong mass concentration in the centre of the disc. Addi¬ 
tional test cases for the numerical method are presented be¬ 
low in the results section. 


4 RESULTS 

4.1 Test case: The pressure-less disc 

Ideally, any newly developed numerical method needs to be 
tested against known analytical solutions. Fortunately, there 
exists a solution in the special case of a fully relativistic 
pressure-free disc (Neugebauer & Meinel 1993, 1995). A se¬ 
quence of dust disc models were run with a varying central 
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Figure 3. The total gravitational mass M versus central redshift 
for the non-rotating disc. The dashed lines refers to the equation 
of state (20) and the solid line to (23). 





Figure 4. The normalized surface density a as a function of 
radius for the non-rotating disc. The labels denote the value of 
the central redshift of the disc. 


redshift. In Fig. 1 the angular velocity fl (solid-line, dia¬ 
monds) and the global quantity Mti (dotted-line, triangles) 
are compared with the analytical solution for different cen¬ 
tral redshift z. The numerical and analytical solutions agree 
very well. Note that the angular velocity is not a given 
fixed input parameter, but is determined self-consistently 
as a function of z (with fixed outer disc radius pa) dur¬ 
ing the iterative solution. The normalized central redshift 
z n = z/(l + z) serves as parameter which measures the 
strength of relativistic effects. In the Newtonian limit z n = 0 
and in the extreme relativistic limit z n = 1. Locally, for the 
metric potentials, the agreement is very good as well. The 
metric potentials displayed in Fig. 2 are those of the line 
element written in Weyl-Lewis-Papapetrou form which has 
been used by Neugebauer & Meinel (1993). They are related 
to v , to, p by 


e 2y = 


2w 2 2 —2u 

e — p to e 


2U 2 —2u 

ae = p cue 

e 2k _ e 2(v+n) 


(36) 

(37) 

(38) 



z/(1+z) 


The potential B is equal to 1 everywhere for the dust 
disc. The results are compared for an intermediate redshift 
z = 2.97. For the general case, with internal pressure, there 
exists no analytical solution for a relativistic disc. In the 
Newtonian limit however, the numerically obtained solutions 
also agree very well with the known analytic Maclaurin disc 
solution (^J, ||) and with the bifurcation diagrams found in 
the Newtonian case (Kley 1996). Although a complete test 
of the numerical approach was not possible, the test cases 
provide very positive indication as to the reliability of the 
results. 


Figure 5. The central value of the potential B versus central 
redshift for the non-rotating disc for two different equations of 
state (see fig. 3). 

4.2 Non-rotating discs 

The non-rotating discs with Q. = 0 are supported purely by 
internal pressure and have, for sequences of a given central 
redshift z, the maximum value of K. Similar to stars, non¬ 
rotating discs have a maximum mass M/pd which increases 
with central redshift. In Fig. 3 the total gravitational mass is 
plotted versus central redshift for the two different equations 
of state. The higher total masses are reached in the case of 
the simple equation (|i^), dashed line, while the polytropic 
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Figure 6. The polytropic constant K versus central redshift for 
the non-rotating disc for two different equations of state (see 
Fig. 3). 

relation (^ ) has a maximum at a finite redshift z « 3.25. 
The stronger relativistic cases have been calculated here us¬ 
ing a higher grid resolution of 256 x 256. However, the last 
part of the curves for z/{ 1 + z) *1 0.9 may nevertheless still 
be unreliable. The problem arises from the strong central 
concentration of the surface density (Fig. 4). Due to this 
problem M reaches its maximum for the first equation of 
state (^) most likely at infinite central redshift. At larger 
and larger redshifts, the density increasingly peakes at the 
origin (refer to the model with a = 17.0). The sequence ter¬ 
minates eventually at a black hole configuration, with all the 
mass located at the origin. The strong increase of internal 
pressure with z leads also to a central drop in the metric 
potential B which is plotted in Fig. 5. For the pressure-free 
dust disc, B = 1 everywhere; however, in the case of non¬ 
zero pressure, equation (^) creates a deviation from this 
simple relation. With increasing z, the metric potential B 
decreases (in all space) and its central value approaches zero 
in the extreme relativistic limit (z = 1). The variation of the 
polytropic constant K with redshift is displayed in Fig. 6. 
With increasing redshift, K drops from the Newtonian limit 
and increases again for stronger relativistic discs. Adding 
rotation (fi > 0) will lower the value of K. In the follow¬ 
ing diagrams K will be normalized to this maximum value 

l/||i;ix ■ 


4.3 Rotating discs 

Rotating discs were studied for the isentropic second equa¬ 
tion of state only (^s|), and the primary results are displayed 
in Fig. 7 in the K — Q diagram. The curves are labeled with 
the corresponding central redshifts of the sequences. Weaker 
relativistic discs (solid line, z = 0.025) follow the Newtonian 
curve (thin dotted line, eq. |rj|) closely up to a finite < f2 c 
and then bifurcate continuously into a ring-like structure at 



Figure 7. Pressure constant K versus angular velocity ft given 
in units of the angular velocity of a dust disc fl c with the same 
central redshift. The solid line refers to a central redshift z = 
0.005, and the dashed dotted line to z = 0.22. The other curves 
are labeled with the corresponding z values. The thin dotted line 
is the Newtonian limit (Maclaurin disc). The thick dotted line 
denotes the mass-shed limit. The values of S7 C and A' max can be 
read off from Fig. 1 and Fig. 6. 



Figure 8. Enlargement of Fig. 7 displaying in detail the bifurca¬ 
tion from the disc to the ring configuration. The dotted line refers 
to the Newtonian disc solution. The solid line refers to the disc 
solution which turns into a ring structure (dashed line) at point 
3. The corresponding surface density distributions at the marked 
locations are displayed in the next figure. 
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Figure 9. Plot of the surface density distribution versus radius 
at the locations marked in Fig. 7. Curve 1 follows closely the 
Newtonian curve 


Cl/Cl c ~ 0.84. This continuous transition from disc to ring 
occurs only for z £ 0.01. The details of the bifurcation pro¬ 
cess will be explained below. For intermediate central red- 
shifts 0.01 < z < 0.22 (dashed-dotted line), the discs and 
rings coexist with no apparent connecting branch between 
them. Stronger relativistic, 2 0.22, (dashed lines) discs 

terminate in a mass shed limit, where at the outer edge of 
the disc (p = pd) gravity is balanced exactly by centrifugal 
forces. From the hydrostatic equation we obtain the criterion 
for reaching the mass-shed limit 


u; p = y, P 


1 V%,p 

2 1 - 1 $ 


= 0 . 


(39) 


The end points of the dashed curves and the thick dotted line 
mark the loci where the edge of the disc reaches this limit. 
In case of the dust disc, all radii are at this limit simulta¬ 
neously, since eq. (^9|) is exactly the hydrostatic equation 
for a pressure-less disc. It should be noted that the angular 
velocity reached in these stronger relativistic cases exceeds 
the value of the pressure-free case ($7 = Cl c ) for the same 
central redshift. This effect is caused by the non-linearity 
of the Einstein-equations. Through the gravitational action 
of the internal pressure, the potential B is then influenced, 
which in turn enters through the rotational velocity v v into 
the hydrostatic equation. 

A more detailed view of the continuous bifurcation pro¬ 
cess is given in Fig. 8, which is an enlargement of the previ¬ 
ous figure. The sequence bifurcates smoothly from the disc 
solution through a sequence of dumb-bell shaped density 
distributions which show a depression of the surface density 
in the middle region of the disc (see Fig. 9). At point 3 the 
density in centre has reached zero, and the ring structures 
begin. The numbers in Fig. 8 refer to the density distribu¬ 
tions displayed in Fig. 9. The solid line in Fig. 9 follows 
approximately the Newtonian relation (ps|). With a growing 
inner ring radius the maximum surface density increases, 
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Figure 10. Detail of the bifurcation process for four different 
central redshifts (dashed line z = 0.0005; solid line z = 0.005, 
the others are given in plot). The Newtonian curve for the disc 
is given by the dotted line. The black squares mark the onset of 
the ring sequences (not plotted for z = 0.015). 


and the ring sequence turns eventually into a line with zero 
radial extent. These extreme ring models could not be ob¬ 
tained numerically, due to the lack of resolution. 

At larger central redshifts, the bifurcation process 
changes its character; the discs and rings become discon¬ 
nected. This is illustrated in Fig. 10, where the bifurcation 
regime in the K — Cl diagram is shown for three different 
values of 2. At a critical value of z, which is given to the 
numerical accuracy by z c ~ 0.01, the connection between 
disc and ring is lost. Approaching z c from below, the curves 
in the K — Cl diagram begin to be double valued until, upon 
reaching the critical point, they disconnect. Above z c , the 
disc sequences do not join smoothly to the ring sequence but 
terminate instead at the mass-shedding limit. The ring se¬ 
quences can still be continued somewhat into the dumb-bell 
regime but there the curve also terminates. Note that the 
apparent difference in the bifurcation curve for the weak 
relativistic limit and the Newtonian limit (see Kley 1996, 
Fig. 1) stems from the fact that here curves of constant 
central redshift (fixed central potential i/(0,0)) were consid¬ 
ered, while in the Newtonian work the curves are plotted 
for constant disc mass. The Newtonian line separates the 
two different regimes. Curves above end in the mass-shed 
limit, while curves below bifurcate continuously into a ring 
sequence. In Fig. 11 it is demonstrated that the disc se¬ 
quences for 2 > 2 C , which terminate at the mass-shedding 
limit, start at that point where the rotational velocity v v at 
the outer rim [p = Pd) coincides with the value of a pressure¬ 
less dust disc having the same central redshift. Above 2 C , the 
rim-velocity lies always below that of the dust disc. 

The total kinetic energy of the disc can be defined for 
rigid rotation by 

T = IciJ. (40) 
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Figure 11. The rotational velocity at the rim (p = p^) of the 
disc as a function of central redshift for the dust disc (dashed 
line) and along the mass-shed limit (solid line). The curves join 
at the point where the disc begin to reach the mass-shed limit at 
z = 0.22. 



Figure 12. The normalized angular velocity versus the ratio of 
kinetic to potential energy for discs with different redshifts. The 
line types and the labeling corresponds to Fig. 7. 

The total proper mass of the system, i.e. the energy of the 
configuration excluding gravitational potential energy and 
kinetic energy, can for the disc be defined as 

M P = / (cr° + (T l )uy/^~gd 3 x. (41) 

■J disc 

The gravitational binding energy is then given by 


W = M p + T - M. (42) 

In Fig. 12 the normalized angular velocity is plotted versus 
the ratio t of the kinetic energy to the gravitational binding 
energy, t = T/\W\. The Newtonian curve (thin dotted line) 
is given by the relation 

Q 2 K= — t(l-2t), (43) 

which terminates at t = 0.5. For stronger relativistic discs 
the maximum possible t becomes smaller and smaller be¬ 
cause the discs reach the mass-shed limit. For the weaker 
relativistic disc the curve for the mass-shed limit approaches 
the Newtonian limit tangentially. 


5 CONCLUSION 

We calculated the structure of uniformly rotating, infinites¬ 
imally thin, relativistic self-gravitating discs with internal 
pressure. The pressure is given by a poly tropic equation of 
state. The polytropic exponent 7 = 3 was used, since in this 
case there exists an exact solution in the Newtonian limit. 
This special choice corresponds to three-dimensional bodies 
of constant densities as well (see Hunter 1972). Thus, we 
may make a comparison of the results found here for flat 
discs with those of rotating, homogeneous relativistic stars. 
As Butterworth & Ipser (1976) showed, sequences of homo¬ 
geneous, rotating relativistic stars usually terminate at the 
mass-shed limit. Close to the Newtonian case however, they 
were not able to follow the sequences to this limit. We sug¬ 
gest that, as weaker relativistic discs bifurcate into a ring, 
weaker relativistic stars of constant density bifurcate into a 
toroid structure as well, and do not end at the mass-shed 
limit. In the Newtonian limit, the structure of these con¬ 
stant density toroids were calculated by Eriguchi & Sugi- 
moto (1981). 

The pressure-less discs possess ergo-regions, where the 
dragging of inertial frames would force observers to ro¬ 
tate. These first appear at a single point within the disc 
at z e = 1.41 and reach the edge of the disc at 2 = 1.89 
(Meinel & Kleinwachter 1993). They also occur in rotat¬ 
ing stars (Butterworth & Ipser 1976). For discs with inter¬ 
nal pressure however, we found no indication of the exis¬ 
tence of ergo-regions. It is to be expected that there exists 
a continuous transition from the zero-pressure K = 0 case 
to discs with non-vanishing pressure (K small). However, 
this connection could not be demonstrated numerically, as 
disc sequences end either at the mass-shed limit or bifurcate 
into rings before reaching the pressure-less limit. One might 
speculate, perhaps, that such a connection can be achieved 
by a disc which consists of two parts: a pressure supported 
central region (0 < p < p p ), surrounded by a dust disc 
( p p < p < pa). The dust disc would refer to p p = 0, and the 
disc sequences to p P = pd- 

In some rotating relativistic stars, there are sequences 
of supra-massive stars (Cook, Shapiro & Teukolsky 1992), 
which are so massive that they exceed the rest mass of a 
non-rotating star, and can only exist for non-zero rotation. 
In the case of flat relativistic discs with non-zero internal 
pressure, we did not find any supra-massive disc sequences. 
Possibly, these negative results are a result of the particular 
equation of state where in the strong relativistic limit the 
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mass is always concentrated in the centre of the disc, as 
shown in Fig. 4 for the non-rotating case. Other relations 
between pressure and surface density, in particular a smaller 
7, could also possibly lead to supra-massive sequences or 
ergo-regions. 

The numerical method developed here is sufficiently 
general to be applied to the study of differentially rotating 
discs, and can easily be extended to three-dimensional rotat¬ 
ing stars as well. This may be the subject of a future further 
investigation. A stability analysis of the computed equilib¬ 
rium configurations lies beyond the scope of the present pa¬ 
per. 


ACKNOWLEDGMENTS 

Useful discussions with Drs R. Meinel and T. Wolf are grate¬ 
fully acknowledged. 


REFERENCES 

Bardeen J. M., Wagoner R. V., 1971, ApJ, 167, 359 
Binney J., Tremaine S., 1987, Galactic Dynamics. Princeton Uni¬ 
versity Press, New Jersey 

Butterworth E. M., Ipser J. R.. 1976, ApJ, 204, 200 
Cook G. B., Shapiro S. L., Teukolsky S. A., 1992, ApJ, 398, 203 
Eriguchi Y., Sugimoto D., 1981, Prog. Theor. Phys., 65, 1870 
Hunter C., Ann. Rev. Ast. & Astrophys., 3, 219 
Neugebauer G., Meinel R., 1993, ApJ, 414, L97 
Neugebauer G., Meinel R., 1995, Phys. Rev. Lett., 75, 3046 
Meinel R., Kleinwachter A., 1995, Einstein Studies 6, Birkhauser 
(Boston), p339. 

Kley W., 1996, MNRAS, 282, 234. 

Pichon C., Lynden-Bell D., 1996, MNRAS, 280, 1007. 

Turakulov Z.Y., 1990, Mod. Phys. Let. A, 5, 725. 


© 0000 RAS, MNRAS 000. 000-000 



